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Abstract 

We have studied the effective response of composites of spherical particles 
each having a dielectric profile which varies along the radius of the particles. 
We developed a first-principles approach to compute the dipole moment of 
the individual spherical particle and hence the effective dielectric response of 
a dilute suspension. The approach has been applied to two model dielectric 
profiles, for which exact solutions are available. Moreover, we used the exact 
results to validate the results from the differential effective dipole approxi- 
mation, recently developed to treat graded spherical particles of an arbitrary 
dielectric profile. Excellent agreement between the two approaches were ob- 
tained. While the focus of this work has been on dielectric responses, the 
approach is equally applicable to analogous systems such as the conductivity 
and elastic problems. 
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I. INTRODUCTION 



Functionally graded materials (FGM) are inhomogeneous materials with spatially vary- 
ing material properties. These materials have received considerable attention as one of 
the advanced inhomogeneous composite materials in various engineering applications since 
first being reported in the 1980s [1]. The change in the composition induces material and 
microstructure gradients, and makes the FGM very different in behavior from the homoge- 
neous materials and conventional composite materials [1,2]. These materials can be tailored 
in their materials properties via the design of the gradients. For mechanical properties, 
the main advantages of a graded material profile range from the improved bonding strength, 
toughness, to wear and corrosion resistance. Other benefits include the reduced residual and 
thermal barrier coatings of high temperature components in gas turbines, the surface hard- 
ening for tribological protection and graded interlayers used in multilayered microelectronic 
and optoelectronic components [1-3]. 

Over the past few years, there have been a number of attempts, both analytical and ex- 
perimental, to study the responses of FGM to mechanical [4-8], thermal [9-11] and electric 
[12,13] loads and for different microstructure in various systems. Nevertheless, the responses 
of composites of FGM inclusion, which is formed when graded inclusions are suspended ran- 
domly in a host medium, should be more useful and interesting. Various different attempts 
have been made to treat the composite materials of homogeneous inclusions [14] as well as 
multi-shell inclusions [15-18]; there exist many methods to study the effective properties of 
composite materials in the literature. However, these established theories for homogeneous 
inclusions cannot be applied. It is thus necessary to develop a new theory to study the 
effective properties of graded composite materials under externally applied fields. In this 
paper, we will develop a first-principles approach for calculating the effective response of 
graded composite materials. 

The paper is organized as follows. In the next section, we present the formalism. We 
will solve the field equations analytically for graded spherical inclusions. In sections III, we 
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obtain the exact solution for a power-law profile, and in section IV, for a linear profile. In 
section V, we use the exact solution to validate a recently proposed approximate theory. 
Discussion and conclusion will be given. 



II. FORMALISM 



In this section, we will focus on the dielectric response as an example. The formalism 
can equally be applied to analogous systems like the conductivity and elastic problems. 
More precisely, for a composite of graded spherical inclusions, we will find its response to 
an externally applied electric field. The formalism involves two major parts: firstly we 
will calculate the local electric field distribution in a graded spherical inclusion and then 
the induced dipole moment of the inclusion. The two-dimensional composite of graded 
cylindrical inclusions has been discussed in [19]. In this work, we will extend the analysis to 
the more realistic case of three-dimensional composites of graded spherical inclusions. We 
will study two model dielectric profiles, namely, the power-law and linear profiles for which 
exact solutions of the local electric field and the dipole moment can be obtained. 

We consider a graded spherical inclusion of radius a suspended in a homogeneous host 

— * 

medium, subjected to a uniform electric field E along the z-axis. We will consider a low 
concentration of suspended inclusions. Thus, the interaction among the inclusions can be 
neglected and the effective dielectric properties of the composite can be obtained from the 
responses of a single inclusion under an effective electric field E. For dielectric response, the 
constitutive relation of a graded spherical inclusion reads 



where 6j(r) is the dielectric profile of the graded spherical inclusion. The relation for the 
host medium is 



where e m is the dielectric constant of the host medium. The Maxwell's equations read 



D = ei(r)Ei. 



(1) 




(2) 
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V-D = 0, (3) 

and 

V x E = 0. (4) 

— * 

To this end, E is the gradient of a scalar potential $: 

£ = -V$. (5) 

Eq. (3) and Eq. (5) can be combined into a partial differential equation: 

V • (e(r)V$) = 0. (6) 

We will normalize the dielectric profile to the dielectric constant of the host medium e m for 
convenience. Without loss of generality, we may also let a — 1. 
In spherical coordinates, the electric potential $ satisfies 

^ ( r e(r) ^J ( Sm ^ (r) ^J + ^in^^ ( e(r WJ = °' ^ 

where e(r) = ej(r)/e m , 6j(r) is the dielectric constant of the inclusion and e m is the dielectric 
constant of the host medium. 

We consider the applied electric field along the z-axis, thus $ is independent of the 
angle if. If we write $ = i?(r)6(6 ) ), after a separation of variables, we obtain the ordinary 
differential equation for the radial function R(r), 

d / 2 dR\ ^ r 2 de(r) dR . _|_ -^\^> q 
dr \ dr J e(r) dr dr 

where n is an integer, e(r) is a dimensionless dielectric constant, while e(r) = 6j(r)/e m in 
the inclusion, and e(r) = 1 in the host medium. 

The potential can be obtained by solving Eq.(8). In the dilute limit, the dipole moment of 

— * — * 

the graded spherical inclusion can be derived. We take the average of the operator D — e m E 
in the whole volume of the composite material, then 
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1 / (D - e m E)dV = D — e m E, (9) 



VJv 

where V is the volume of the whole composite material, A denotes the average of an operator 
A in the composite material. The integrand vanishes in the host medium, and thus Eq.(9) 
becomes 

yJ Q (D- e m E)dV — D — e m E, (10) 

where Qi is the region occupied by the inclusion. 

Now, we can define the effective constitutive relation for the composite material: 

D = e e E, (11) 

where e e is the effective dielectric constant of the composite material. Thus 

\- I (q - e m )EdV = (e c - e m )E. (12) 

Eq.(12) gives the polarization of the composite and it can be used to calculate the effective 
dielectric properties of the composite material at a low particle concentration. 

III. EXACT SOLUTION FOR A POWER-LAW PROFILE 

The dielectric constant in the particle is taken as a power-law function of the radius. In 
this case, e(r) = c^r^, with k > where < r < 1. Then the radial equation becomes 

d 2 R k + 2dR n(n + l) 
dr 2 r dr r 2 

As this is a homogeneous equation, it admits a power-law solution, 

R(r) = r s . (14) 
Substituting it into Eq.(13), we obtain the equation 

s 2 + s(k + l) -n(n + l) = 0. (15) 



The solution to this equation is 

, 1 



(k + 1) ± yj (k + l) 2 + 4n(n + 1) 
The potentials in the inclusion and the host medium are given, respectively, by 

oo 

$,(r, 9) = A + E [A„r s +^ + B n r s -^] P„(cos#), 

n=0 

oo 

$ m (r, 9) = C + J2 [C n r n + A>r- (n+1) ] P n (cos 9). 



n=0 



As r — > oo, the potential is given by the far field 

$ m = -E r cos 9, 

while r — > 0, $j attains a finite value. Thus 

$i(r,0) = A 1 r s + (1) P l {cos9), 
$m(r,6) = -E r cos 9 + ^P 1 (cos 9). 

Meanwhile, the potential functions should satisfy the boundary conditions, as follow: 

$i(r,0) | r= i = $ m (r,0) \ r =i , 
d$i{r,6) . 



,a$ ro (r,0) | 



dr 



r=l 



We obtain: 



-E Q + D X = A X , -E -2D 1 = c k A 1 s k + (l). 



It is not difficult to solve these equations and obtain the coefficients, 



-_E , Z^! — ; ^ E . 



sc k + 2 



sc fe + 2 



where s = s%(l) = § y (A; + l) 2 + 8 - (fc + 1) . Then the potentials are given by 

$ m (r, 9) = -E r cos 9 + ^ cos 0, 
$i(r, 0) = Air s cos0. 



The electric field E{ along the z-axis in the inclusion can be derived from the potential 



E iz = -A x r 



s-l 



's - l)cos 2 ^ + 1 



(24) 



In the dilute limit, the polarization can be calculated from the following equation: 



ATre m E b = (e^r) - e m )E iz (r, 6)dV, 



(25) 



where b is the dipole factor, which measures the degree of polarization of the graded spherical 
inclusion against the host medium, and Qi is the region occupied by the inclusion. After 
performing the integral, we obtain 

s + 2 / c k 1 



b = 



scfc + 2 \k + s + 2 s + 2 
When k — 0, e(r) = Ch is a constant, then 

s = 1 



(26) 



(27) 



and 



c fc + 2 

This is just the dipole factor in the homogeneous case. 



(28) 



IV. LINEAR PROFILE WITH A SMALL SLOPE 

In this section, we consider e(r) as a linear function of the radius. This is a simple 
graded material that can be made easily. Moreover, other materials can be simulated using 
piecewise linear functions. For the inclusion, 

e(r) = d + cr. (29) 

According to Eq.(8), the radial function for the inclusion satisfies the following equation: 
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SR , (2 , 1 ^ n(»+l) fl = a (3Q) 



df 2 Vr f + 1/ df r 2 
where r = £f, and £ = d/c. We consider the case with a small slope, that is, |£| = |d/c| > 1, 
we can express the solution in a power series. The power series solution is: 



fn(r) = £ C%r k+p . (31) 

fc=0 

Then the radial equation [Eq.(30)] becomes 

oo 

£ C£ [(/c + p)(k + p - 1) + 2(k + p) - n(n + 1)] r fc+ ^ 2 

fc=0 

oo 

+ E Cfc [(* + P)(A: + p - 1) + 3(A; + p) - n(n + 1)] f^ 1 = 0. (32) 

fc=0 

The coefficient of each term should vanish, and the lowest term satisfies: 

[p(p - 1) + 2p - n{n + 1)] CI = 0. (33) 

This equation is for the characteristic equation of the differential equation [Eq.(30)]. Solving 
it, we obtain 

p = n or — (n + 1). (34) 

Similarly, the recursion relation can also be found from Eq.(32), 

_ (k + n)(k + n + 2) - n(n + 1) 
fc+1 (A; + n + l)(A; + n + 2)-n(n + l) fe ' 1 J 

This power series solution is absolutely convergent in the region < f < 1. Thus, we can 
use this solution to study the response of composite material to the external field, as long 

as \d/c\ > 1. 

Consider the properties of the potential at infinity and at the origin, the potentials are 
given by 

$ m (r, 9) = -E r cos 9 + ^ cos 9, 

oo 

$i(r,e) = A 1 ^C 1 k (-r) k+1 cos9. (36) 
k=o d 



Moreover, from the boundary conditions, we obtain the coefficients A\ and D±, 

-3E (d + c)v 2 -vj_ 

Al ~ (d + c)v 2 + 2v l > Dl ~ (d + c)v 2 + 2v 1 E °- (37) 



where 

fc+l oo / r \ k+1 

) , ^2 = E < ' 

k=0 k=0 

Then the z-component of the electric field can also be calculated 

fc+i 

k=0 

Thus we obtain the dipole factor 



"■ = £3 5 ) , "* = £<*(* + !) 5 ) . (38) 



00 /r\ fe+1 

= -Ai E (^J r fc (A;cos 2 ^ + l). (39) 



6 = + 9„ [(& " 1)Vl + • (40) 

{a + cj-u 2 + 2fi 



where 



«. = £Ms) *+? < 41 > 

When c = 0, e(r) is a constant, Eq.(40) reduces to the homogeneous case. Thus, when 
c — > 0, the dipole factor is given by 

6 = d^l = (42) 
This result is the same with that in the homogeneous case. 

V. COMPARISON WITH DIFFERENTIAL EFFECTIVE DIPOLE 

APPROXIMATION 

In this section, we compare the exact results with the differential effective dipole approx- 
imation (DEDA) [21], recently developed to treat arbitrary graded profiles. We start with 
the differential equation for the dipole factor [21]: 

f r = -3^^) K 1 + 2b ^ ~ (! - 6 ) e ( r )H(l + ^)e m + 2(1 - 6)e(r)], (43) 
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where e m = 1. Thus the dipole factor of a graded spherical inclusion can be calculated 
by solving the above differential equation with a given graded profile e(r). The nonlinear 
first-order differential equation can be integrated, if we are given the graded profile e(r) and 
the initial condition b(r = 0). 

We have evaluated DEDA for for two model graded profiles: (a) power-law profile e(r) = 
Cr n , and (b) linear profile e(r) = A + Br. Note that we have changed notations slightly to 
agree with those of Ref. [21]. The numerical integration has been done by the fourth-order 
Runge-Kutta algorithm with a step size Sr = 0.01. In Fig. 1(a), we plot the dipole factor b 
versus C for various index n > 0. It is clear that b increases monotonically as the dielectric 
contrast C increases, while it decreases with the index n. It is attributed to the fact that the 
average dielectric constant decreases as n increases. Similarly, in Fig. 1(b) we plot b versus 
A for various slope B. We obtained similar behavior as in Fig. 1(a). 

It is instructive to compare the exact results with the DEDA results. In Fig.l, we 
compared the exact results with the DEDA results for a spherical inclusion with a power- 
law profile and a linear profile. The agreement is excellent. Note that exact results are 
unavailable for A < B in the linear profile. Thus, we would say the DEDA is a very good 
approximation for graded spherical inclusions. 

VI. DISCUSSION AND CONCLUSION 

Here a few comments are in order. In this work, we have developed a first-principles 
approach to compute the dipole moment of the individual spherical inclusion and hence 
the effective dielectric response of a dilute suspension. The approach has been applied to 
two model dielectric profiles, for which exact solutions are available. Moreover, we used 
the exact results to validate the results from the differential effective dipole approximation, 
which is valid for graded spherical inclusions of an arbitrary dielectric profile. Excellent 
agreement between the two approaches were obtained. Note that an exact solution is very 
few in composite research and to have one yields much insight. It is thus worth spending 
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time on finding one. To this end, it is also instructive to obtain the analytic result for the 
piecewise linear profiles [19]. 

It is instructive to extend the present approach to nonlinear graded composites. The 
introduction of a graded dielectric profile in nonlinear composites provides an extra dimen- 
sion to tune the enhanced nonlinear response. For instance, a graded interfacial layer on the 
spherical inclusions may help to control the local field fluctuations, and hence the nonlinear 
response. The perturbation approach [22] as well as the variational approach [23] are just 
suitable for this extension. 

Our approach can be applied to biological cells, as the interior of biological cells must 
be inhomogeneous in nature and can be treated as a graded material. To this end, Freyria 
et. al. [24] observed a graded cell response when they studied the cell-implant interactions 
experimentally. Thus the complex dielectric function can be modeled to vary continuouly 
along the radius of the cells, namely, the conductivity can change rapidly near the boundary 
of cells and a power-law profile prevails, while the dielectric function may vary only slightly 
and thus a linear profile suffices. Work is in progress along this direction [25]. 
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FIGURES 

FIG. 1. The dipole factor b plotted for two model profiles: (a) versus C for various index n in 
the power-law profile, and (b) versus A for various slope B in the linear profile. 
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